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Abstract 

■ We analyze a single-spin measurement using a transient process in 

magnetic force microscopy (MFM) which could increase the maxi- 
Q_i| mum operating temperature by a factor of Q (the quality factor of 

the cantilever) in comparison with the static Stern-Gerlach effect. We 
obtain an exact solution of the master equation, which confirms this 
^ | result. We also discuss the conditions required to create a macroscopic 

Schrodinger cat state in the cantilever. 



1 Introduction 

The routine magnetic force microscopy (MFM) seems to be of a little use 
for a single spin detection in solids. Even for unrealistically small distances 
between the ferromagnetic particle on the cantilever tip and the spin, the 
maximum temperature for a single-spin measurement comes to a millikelvin 
region [1]. In this paper, we show that using a transient process, one can 
increase the maximum temperature of a single-spin measurement by a factor 
of Q (the quality factor of the cantilever). Alternatively, one can increase 
the distance between the ferromagnetic particle and the cantilever tip if one 
is willing to work at millikelvin temperatures. 
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In section 2 we explain the basic idea of our work. In sections 3-5 we obtain 
an exact solution of the master equations, which confirms our idea. In section 
6 we derive the conditions for creating a macroscopic Schrodinger cat state 
(MSCS) in MFM. In Conclusion we summarize our results. 



2 Transient Process in MFM 

In this section we describe the basic idea of our work. Suppose that a ferro- 
magnetic particle on the cantilever tip interacts with a single spin on the solid 
surface. (See Fig. 1.) The equilibrium position of the cantilever tip depends 
on the spin direction. The distance between two possible equilibrium posi- 
tions, corresponding to two spin stationary states, is given by 2F/ k c , where 
k c is the cantilever spring constant, F is the magneto-static force between 
the ferromagnetic particle and the spin. In order to measure the state of 
the spin, this distance must be greater than twice the uncertainty due to the 
thermodynamical noise of the cantilever position. This uncertainty can be 
estimated as (A; B T//c c ) 1//2 , where fee is the Boltzmann's constant and T is the 
temperature. Thus, the condition for a single spin measurement (the static 
Stern-Gerlach effect) is [1] 

T <T max =F 2 /k B k c . (1) 

Now, assume that we quickly change the stationary state of the spin and 
consider the transient cantilever vibrations after this change. The amplitude 
of the cantilever vibrations at time t <C Q /uj c (oo c is the cantilever frequency 
and Q/ui c is the time constant of the cantilever) is 4F/k c . In addition, as- 
sume that we detect the position and momentum of a point on the cantilever 
tip with an accuracy that satisfies the quantum limit (5P Z )(SZ) ~ h/2 (the 
cantilever oscillates along the z-axis.) To find out the state of the spin we 
are going to compare the observed trajectory of the cantilever tip with the 
theoretical prediction. The theory predicts the cantilever trajectory within 
an uncertainty due to the thermal noise. However, the thermal uncertainty of 
the cantilever position at time t <C Q/oj c increases as [t(uj c /Q)(k B T/k c )} 1 ^ 2 . 
We can obtain this expression assuming an initial (t <C Q/oo c ) thermal dif- 
fusion with two common properties: 1) The dispersion (SZ) 2 is proportional 
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to time, t, and 2) The uncertainty of the cantilever position equals its ther- 
modynamical value if we formally put t = Q/u c . 

At the time t = ti/uj c (half of the first period), the distance between two 
possible cantilever positions takes its maximum possible value 4F/ k c . At the 
same time, the thermal uncertainty of the position predicted by the theory 
is still much smaller than its thermodynamical value. Now, the condition for 
a single spin measurement is 

T < Tmax = ^kj c - (2) 

One can see that the maximum temperature for a single spin measurement 
increases by a factor of Q, compared with the static Stern-Gerlach effect. 

In the next three sections, we will confirm this estimation by direct solu- 
tion of the master equation. 



3 The Hamiltonian and the Master Equation 

We assume that the cantilever tip with an attached ferromagnetic particle 
can oscillate along the z-axis. (See Fig. 1.) A single paramagnetic atom with 
spin 1/2 is placed near the cantilever tip. The dimensionless Hamiltonian of 
the cantilever tip interacting with a single spin is 



H =-(?, + ?)- 2nzS,. 



(3) 



We introduced the following notation (below we omit hats for operators): 

dB 



Z/Z q) p z = P z /P q , r] = 9f^E 



dz 



(4) 



where Z q , P q and F q are the "quants" of the coordinate, momentum and 
force acting on the cantilever, 



TlUJr 



1/2 



H/Zg, 



F q = k c Z q . 



(5) 
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The variables Z and P z are the "dimensional" coordinate and momentum of 
the cantilever tip, k c is the cantilever "spring constant" , oj c is its frequency, 
g is the "g-factor" of the spin (below we use g — 2), dB z /dz is the magnetic 
field gradient produced by the ferromagnetic particle at the spin location 
when the cantilever is in the equilibrium position with no spin [z — 0). Note 
that the cantilever interacts with the z-component of the spin, which is an 
integral of motion in our system. In the Hamiltonian (3) we omitted the 
term {gix B BQ/huj c )S z , where B is the magnetic field on the spin when the 
cantilever is at the origin {z — 0). This term may be eliminated "physically" 
(by application of a uniform external field of magnitude B$ in the negative z- 
direction) or "mathematically" (by transferring to the system of coordinates 
rotating with the frequency g/i B B /h). 

The master equation describes the evolution of the density matrix of the 
system interacting with the environment (see for example [2-5]). We are tak- 
ing into account the interaction of the cantilever with its environment, and 
ignore the direct interaction between the spin and the environment, assum- 
ing that the spin relaxation and decoherence times are large enough. The 
effect of the environment depends on its "spectral density", i.e. the den- 
sity of environmental oscillators at a given frequency. Probably, the simplest 
model of the environment is the "ohmic" model, where the spectral density 
is proportional to the frequency ui for u < Q, where Q is the cutoff frequency 
for the environment. A master equation for the non-ohmic environment has 
been derived in [5]. For the ohmic model, the simplest master equation has 
been obtained in [2]. This equation is valid in the "high temperature limit" 
k B T » HQ. The master equation derived in [3] is valid for arbitrary temper- 
ature. As pointed out in [5], both equations [2] and [3] fail at times shorter 
than or close to h/k B T. 

We are going to consider the "gedanken experiment" discussed in [1]. Sup- 
pose that initially (t = 0) the spin is in a superposition of the two states with 
the z-projection of the spin S z — 1/2 and S z = —1/2. These two states of the 
spin correspond to two different equilibrium positions of the cantilever tip. 
Thus, the cantilever (without decoherence) would transform into a MSCS: 
two simultaneous equilibrium positions. Certainly, decoherence will destroy 
this state. The master equation describes both the appearance of the MSCS 
and its destruction due to decoherence. 
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Following [1], we consider the ultra-thin cantilever reported on [6]. It has 
the spring constant k c = 6.5 x 1CT 6 N/m, the frequency, uj c /2tt = 1.7 kHz, 
and the quality factor, Q = 6700. The ferromagnetic particle on the can- 
tilever tip is taken to be a sphere of radius R = 15 nm at a distance 5 nm 
from the paramagnetic atom. (Below we consider conditions for increasing 
the distance between the cantilever and spin.) For this case, the static dis- 
placement of the cantilever tip due to its interaction with the single spin 
exceeds the thermal vibrations of the cantilever for temperatures 

T<T_=<^,1.7rf. (6) 

In our gedanken experiment for the temperature T » huj c / k B m 8 x 10~ 8 K, 
we can use the simplest high temperature limit in the "ohmic model". 

The master equation in the high temperature limit can be written in the 
form [2] 



dps 



~5F 



| (d zz - d z , z ,) - J(^ 2 - z") - \{3{z - z')(d z - d z ,) - D(3(z - z>) 



l\2 



—2ir)(z's' — zs) 



(7) 

Here, s and s' take values ±1/2 (we use s instead of S z ), r = u c t, /3 — 1/Q 
and D = k B T/hu c . Using new coordinates 



r = z-z', R = -(z + z'). 

2 



(8) 



Eq. (7) can be written as: 

dp a , s >(R, ryr) 



idji r — iRr — f3rd r — Dj3r 2 — irj 



(2R-r)s f - (2R + r)s 



p s ,s>(R, r,r). 

(9) 
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Performing a Fourier transformation of this equation with respect to the 
variable "R" , one obtains, after re-arrangements, 



dr 
where 



dr 



+ 



r+2r](s'-s) 



dps 



dk 



-D/3r 2 +ir]r(s'+s) 



Ps,s>, (10) 



/+oo 
e tkR p s>s ,(R,r,r) dR. 
-oo 

We can study separately the spin diagonal case (s = s') and the off-diagonal 
case (s 7^ s'). For s' = s (up-up or down-down spins), we have the following 
equation: 



dp s ,s . fa ix d Ps,s . dp s 

+ \pr — k) - + r- 



dr dr dk 

and for s' ^ s (up-down or down- up spins): 



-Df3r 2 + 2irjrs ) p 



dps 



-Dfir 2 ps^s 



(11) 



(12) 



dr dr dk 

We will derive the exact solution of the master equation (7) for the case when 
the spin is "prepared" initially in the superposition of two states with s = 1/2 
and s = —1/2, while the cantilever tip is in the quasiclassical coherent state 



(tt)V4 



exp[ip z - (z - z ) 2 /2] 



(13) 



where the amplitudes a and b correspond to the values of s = 1/2 and 
s = —1/2 respectively. The corresponding density matrix can be written as 

p S siz, z', 0) = 4>(z, s, 0) ® ^(z', s', 0). (14) 

Note that we consider an ensemble of spin-cantilever systems with the same 
initial state. This implies that the experimenter can detect the position and 
momentum of a point on the cantilever tip with quantum limit accuracy 
(Sp z ) 2 (Sz) 2 = 1/4. (In our gedanken experiment, this corresponds to an 
uncertainty of 300 fm for position and 300 nm/s for velocity.) Based on the 
master equation, we can predict the average position of the cantilever tip for 
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its given initial state, depending on the spin state. If the double uncertainty 
of the position is smaller than the separation between two possible average 
positions, the cantilever tip will measure the state of the spin. 



After Fourier transformation, the "cantilever part" of the density matrix 
is represented by 



p S;S ' (k, r, 0) oc exp ip r + ikz — r 2 / 4 — k 2 /A 



(15) 



4 Solution for spin diagonal matrix elements 



The equations for the characteristics of Eq. (11) are 

dr dk dp SyS i 



dr 



or, explicitly 



Pr-k 



dr 
~dr 

dk 



-Df3r 2 + 2ir/sr ) p. 



fir — k, 



r, 



dp s 



dr 



-D(3r 2 + 2ii]sr ) p 



From the first two equations in (17), one obtains 

d 2 k dk 

d^- (3 d-T +k = > 

which has the following general solution 

k = e PT/2 (ci cos 9t + c 2 sin Or^j , 



(16) 



(17) 



(18) 



where = y 1 — (Here we are considering the case (3 < 2 so 6 is a real 
number. The case (3 > 2 can also be solved analytically.) Using the second 
equation in (17) one obtains: 



oPr/2 



cos 9t — 9 sin Or j c x + (j- sin Or + 9 cos 9r j c 2 



(19) 



Inverting Eqs. (18) and (19) as a function of C\ and c 2 one obtain the char- 
acteristic curves: 

c 1 = e-^ 2 (q 1 k + q 2 r), (20) 

and 

c 2 = e-^^( Pl k + p 2 r), (21) 
where the time dependent constants q 1 , q 2 , P\ and p 2 have been defined as 

qi — ^ (J^ sin 9t + 6 cos 9r^j , 

<?2 = -g sin 6>r, 

Pi — ^ cos9t + 9sm9r^, 

P2 = \ cos 9t. 



Substituting (19) into the third equation of (17) and integrating in time, one 
obtains: 



(22) 



p s , s (k, r, t) oc Q(a, c 2 ) exp i2r]s(c 1 g 1 + c 2 g 2 ) - D(3(c\fi + 2cic 2 / 3 + c 2 2 f 2 ) 



where the functions //s and g[s are defined as 

f i ( r ) = (j3 + + p cos 2#r - 2# sin 29r 

f 2 ( r ) = (/? + M!) - /J cos 2#r + 2# sin 29t 



(23) 



/a(r) = 



29 cos29t + /3sin20r 



(24) 



9i (t) = e^ 2 cos 0t, 
g 2 ( T ) = e^sinflr. 



The arbitrary function A which depends on the characteristics is determined 
by the initial density matrix p StS (k(0), r(0), 0), 



Mci, c 2 ) = Ps,s (ci, |/3ci + 9c 2 , O^j exp [-2irjs(cig w + c 2 g 20 )} 
Dp(cjf 10 + 2cic 2 / 3 o + clf 20 )\ , 



x exp 



(25) 



where fio = fi(0) and g^ = gi(0). From the initial density matrix (Eq. (15)), 
we obtain 



p S:S > (k, r, 0) oc exp {i [(|p /? + *o) Ci + p #c 2 } 

xexp{-[(S + |)c? + ?c lC2 + ^]}. 
Substituting (25) and (26) into (23) one obtains: 



(26) 



p S)S (fc,r,r) ocexp^i 



x exp 



-p f3 + z + 2r ] sG 1 ) ci + ( Po 6 + 2r]sG 2 ) c 2 



'P 2 A 2 P0 8 2 2 

,i6 + 4 Cl + T ClC2 + T C2 



x exp {-D(3{F lC \ + 2c lC2 F 3 + F 2 c|)}, 

where Fi and are defined as 

Fi(r) = fi(r) - f i0 , Gi(r) = g^r) - g i0 . 



(27) 



Substituting in (27) the values of characteristics as a function of k and r 
(Eqs.(20) and (21)), one obtains: 



p aa {k, r, r) oc exp -r 2 d + irC 2 + (iB 2 - rSi) - a;k 2 , (28) 



where 



a 2 = e-» 



'P 1 
16 4 



I 2 
4 



(29) 
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B l = e-^|(S + \) 2q ig2 + f(q lP2 + q 2Pl ) + f 2 Pl p 2 
+2Df3[F 1 q 1 q 2 + (q ± p 2 + q 2 Pi)F 3 + F 2Pl p 2 } \, 



(30) 



B 2 (s) = e~^ 2 



-poP + z + 2r 1 sG 1 ) qi + ( Po e + 2r]sG 2 ) Pl 



(31) 



Ci 



-/3t 



(S + I) «2 + ^ + U 



+D(3{F iq 2 2 + 2g 2 p 2 F 3 + F 2 pj) 



(32) 



C 2 (s) = e"^ 2 



PoP + z + 2r]sG 1 I q 2 + (p 9 + 2r]sG 2 )p 2 



(33) 



In Eqs. (31) and (33) the explicit dependence on s is presented. Performing 
the inverse Fourier transform one obtains 



Pi/2,i/2(R,r,r) 



( exp r(\ ■ irC 2 (l/2) 
* x exp [(-rB 1 + iB 2 (l/2) - ifl) 2 /^], 



|6 



P-i/ 2 ,-i/2(R,r,T) = -P— exp -^ + ^(-1/2) 

v 71 " °"* 



x exp [(-rSi + iB 2 (-l/2) - ii?) 2 /4a 2 ]. 



Eqs. (34) represent two squeezed Gaussians with modulus 



(34) 



\Pi/2,i/2(R, r, t) 



exp 



xexp [-(5 2 (l/2)-i2) 2 /4(7*], 



\b\'' 



(35) 



|p-i A -i/ 2 (it:,r,T)| = ^=^ exp [-r 2 (d - S 2 ^ 2 ) 

*xexp [-(B 2 (-l/2)-RY/Aa% 
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Fig. 2 shows schematically two peaks (seen from the top as ellipses) corre- 
sponding to the two matrix elements lpi/2,1/2 1 and |/?_i/ 2> — 1/2 j ■ We denote the 

centers of the ellipses, which lie on the diagonal z = z', by M ++ and M , the 

semi-major axis by a d , the semi-minor axis by a' d and the distance between 
the centers by A d . The ptfi% 1/2 is located at M ++ = (r = 0, R = B 2 (l/2)) or 
z = z' = £2(1/2), while the p™?J 2 _ 1/2 is at M__ = ( r = 0,R = B 2 (-l/2)) 
or (z — z' — B 2 (— 1/2)). The distance A d is given by 

A d = £3(1/2) -B 2 (-l/2). (36) 

From Eqs. (35), we obtain = \[2 a*, and 

2 - = (37) 

For a single spin measurement, the two peaks corresponding to p r ^fi/2 anc ^ 
P-1/2 -1/2 m ust be well separated. It follows that the condition A d > 2<j d 
must be satisfied. 

First, we consider the case (3t ^> 1 or r ^> Q/u c , where Q/u c is the time 
constant for the cantilever. In this case, we obtain two equilibrium positions 
for the cantilever, when the transient process is over. We have A^ = 2i] 
and o& = \f~D. The value a d = V^D is the thermodynamical uncertainty 
in the cantilever position caused by the thermal noise. The two equilib- 
rium positions can be distinguished if 77 > VD or T < F 2 /k B k c , where 
F = p B \dB z /dz\ is the magneto-static force between the ferromagnetic par- 
ticle and the paramagnetic atom. The last expression exactly coincides with 
formula (1). 

Next, we consider the initial transient process after the instant (t = 0) at 
which the paramagnetic spin has been transferred into the superpositional 
state. For f3r <C 1, we have 

A d = 4r7sin 2 ^, a d = [l/2 + Dfir - Dficosr sin 3 r] 1/2 . (38) 

This expression for A d describes the oscillating distance between the two 
peaks. It corresponds to initial vibrations of two classical oscillators near 
their equilibrium positions z — rj and z = —rj. The distance between them 
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is given by A<j. (For our gedanken experiment the maximum value of A^ 
is 0.24 nm.) The formula for ad contains three terms. The first term, 1/2, 
corresponds to the quantum dispersion of the initial wave function. The sec- 
ond term, D(3t, describes the initial diffusion of an ensemble of oscillators. 
Formally, setting r ~ 1//3, we can estimate the final dispersion = \/~D, 
which corresponds to thermodynamical vibrations of the cantilever tip. The 
third term describes insignificant oscillations with small amplitude, Df3. 

Note that the condition for distinguishing two cantilever positions at the 
beginning of the transient process is much less restrictive than the corre- 
sponding condition for the equilibrium positions at /3t >> 1. Indeed, after 
the first half-period (r = 7r), we have A^ = An and ad = (1/2 + nDf3) 1 ^ 2 . 
Taking into account that /5 1, the condition for distinguishing two posi- 
tions, n > (1/2 + irDft) 1 / 2 , is much easier than rj > \f~D. In our gedanken 
experiment the condition for distinguishing the two positions for the transient 
process is 

40F 2 

T <T max = -?— = UK, (39) 

B C 

compared with T < T max = 1.7 mK for the static Stern-Gerlach effect. This 
estimate seems to be too optimistic. It is connected with the very small 
distance (5 nm) between the ferromagnetic particle and the paramagnetic 
atom. If we increase this distance to 50 nm, the temperature T max drops 
from 14 K to 1.1 mK. Note that expression (39) coincides exactly with our 
preliminary estimates (2). 

The condition A d > 2a d is satisfied for the first time at 

r = T w 2 1 /4/ v ^. (40) 

This expression is valid if n >> 1 and n >> (D(3) 2 /y/8. For our gedanken 
experiment we have n = 144, D = 1.25 x 10 7 T (T is the temperature in 
Kelvin), f3 — 1.5 x 10~ 4 , and r = 0.1. Thus, the above conditions are both 
satisfied. The value of to = uj c t is approximately 9.3 /is. 
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5 Solution for off-diagonal spin matrix ele- 
ments 



The equations for the characteristics are now given by 

dr dk dp s - s 



or 



dr = 



(3r — k r — Arjs —D(3r 2 p s 



dr 



M = r _ 4r?S; 

%i = -D^ p at _ a . 
The solutions of the first two equations of (42) are 



(41) 



(42) 



k = e^ T l 2 ( a cos 9t + c 2 sin Or ) + Af3rjs 



(f cos 6t -6 sin 0r) c x + (§ sin 9t + 9 cos 6>r) c 2 



+ Arjs. 



(43) 



Following the same steps as above we obtain for the Fourier transform: 



Pi/2,-i/2(k,r,r) cx A(d,c 2 )exp^-DP fic\ + 2c l c 2 f 3 + f 2 cl j 



x exp j-£>/3 



^(s-iCi + # 2 c 2 ) + 4?7 2 r 



(44) 



where we fix s = 1/2 (changing sign of s corresponds to a change of sign 
of r), see Eq.(42), therefore the case s = —1/2 can be easily obtained). The 
functions fi and ^ have been defined as above, and c\ and c 2 are new char- 
acteristic curves given by 



Ci =e P T/2 (qik + q 2 r + r/q 3 ), 
(pik + p 2 r + f]p 3 ). 



c 2 



-fir fi 



(45) 
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Here, q±, q2,Pi,P2 are defined by Eqs. (22) and ^3,^3 are given by 



?3 



P3 =1 



(3 § sin 9t + 9 cos 9t + sin 9r 



P ^| cos 9r — 9 sin 9r^j — cos 9r 



(46) 



With the same initial condition Eq. (15), we can determine the function 

^4(ci,c 2 ) and obtain 



P1/2 -i/2(k, r, t) oc pi/2 _i/2 ^ci + 2/3?7, |/3ci + #c 2 + 2?y, 0^) 



x exp <^ -DP 



F lC \ + 2cic 2 F 3 + F 2 c\ + 4r7(GiCi + G 2 c 2 ) + 4t? 2 t 

'(47) 

where Fj(r) and Gj(r) are defined as above. By substituting the initial 
condition (15) we have 



Pi/2,-1/2 (A;, r,r) oc exp <j i 



-Pop + z )d + p 9c 2 + 2r](p + z p) 



x exp 

x exp < —^4:1] 
x exp <^ -DP 



2 80 Pc\? 
16 + 4 Cl + T ClC2 + ^ 



+ v 2 (1 + /? 2 )}} 

F lC 2 + 2ciC 2 F 3 + F 2 c\ + 4r7s(GiCi + G 2 c 2 ) + ?? 2 t 



}• 

(48) 



which can be written as 



Pi/2,-i/2(fc, r) oc exp [-r 2 C 12 - njCn - r] 2 C w + irC 2 \ + ir)C 20 
x exp [(iB 20 - rB n - rjB 10 ) k - a 2 k 2 }, 

where a* is given by (29) and 



(49) 
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c 



12 



,-0T 
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-/3t 



+e 



-/3r 



^ + | + D/JFxj 2g 2 g 3 + + 2L>/3F 3 ) (g 2 p 3 + p 2 <7 3 ) 
€ + D(3F 2 \ 2p 2 p 3 



+4e -/3r/2 



2 )Vl 



Cio =e"^ 



§ + J + D/JFx) g 3 2 + f ^ + 2L>/3F 3 ) g 3 p 3 +(£ + D(3F 2 \ p\ 



+4e~ /3r/2 



+ DPd) q 3 + (^ + D{3G 2 ) p 3 



C21 
C20 



= e 



-0r/2 



^ + Zoj <?2 + PO0P2 



+ 2{p + z P) , 



B 
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-0T 



+e 



(C + 1 + 2 92<7i + + 2£>/3F 3 ) (g lP2 + q 2Pl ) 

£ + D[3F 2 ) 2p 2Pl 



B 
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-fir 



+e 



T 



+ D(3F 2 2p 3Pl 



+4e -/3r/2 



^ + DPG 1 ) qi +(i + D[3G 



2 Pi 



B 20 =e^/ 2 



(]}PoP + zo) q\ +Podpi 



Performing the inverse Fourier transform we obtain 



(50) 
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Pi/2,-i/2(-R, r,r) 



P-i/2,i/2(R,r,r) = 



ab* 



x exp 
a*b 

\fncr* 
x exp 



cxp 



-r 2 C 12 - rrjCu - r/ 2 Ci + irC 2 i + irjQ 



20 



-rB u - r)B w + iB 20 - iRf /4o* 



exp 



-r 2 C 12 + rrjC n - rj 2 C 10 + irC 21 - irjQ 



20 



(-rfln + r]B 10 + iB 20 - iR) 2 j\a\ 



Taking the modulus one obtains 



(51) 



\Pl/2,-l/2(R, T,T) 

|p-i/2,i/2(-R, r,r) 



where 



e ^ e -(r + r 0?7 ) 2 /2a 2 e -(S 20 - i2) 2 /4o* 

\/7rcr* 

M_ e C, 2 e -(r - r r/) 2 /2^ 2 e -(5 20 - i?) 2 /4<7 2 
\/7ro"* 



(52) 



2aj 

4^C 12 -B^' 



ii 



BttB 



11-DlO 



(53) 



t - 5 io 



7 



10 



2a z 



The maxima are located at [R = B 20 ,r = —rjro) for |pi/ 2 ,-i/ 2 |) and at 
(R = B 2Q ,r = r/ro) for |p-i/ 2 ,i/ 2 |. In (z,z') coordinates this corresponds 
to M + _ = (z = B 20 - 7770/2, z' = .620 + r/r /2) for |pi/ 2 _i/ 2 | and M_ + = 
[z = B 20 + 7770/2, 2/ = _B 20 — r/r /2) for |p-i/ 2 ,+i/ 2 |, so that the distance be- 
tween them is given by A n ^ = \f2ri\ro\. Next, we consider the quadratic form 
(r ± r r/)/2cr 2 + (B 20 — R) 2 /<i<j 2 in the (z,z f ) plane. Straightforward calcu- 
lations show that this is an ellipse whose semi-axes are respectively given by 
a (across the diagonal) and 2\f2o* (along the diagonal). The centers of the 
peaks M + _, M_ + are symmetric with respect to the diagonal line z = z' . 



The most remarkable difference compared with the diagonal case, is the 
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presence of irreversible decoherence. Indeed, the heights of the peaks are ex- 
ponentially reduced in time by the damping factor ~ exp(— 4i] 2 Dj3r). This, 
in turn, defines a characteristic time scale of decoherence: r d = 1/At] 2 D(3. 
This formula exactly coincides with the expression derived in [1], based on 
a semi-qualitative analysis. At time r = r , when two diagonal peaks are 
clearly separated, the damping factor is Ai] 2 Dj3TQ. We expect to observe the 
coherence between the two peaks (MSCS) if this factor is not much more 
than one unit. Thus, using the expression for r from Eq. (40), we can 
estimate the condition for the quantum coherence as Df3rfl 2 < 1, or 

max = Y B — W 2 — ' ( ' 

For our gedanken experiment the value of T max is approximately 3 x 10~ 7 K. 

Now, we will check the validity of our estimate for the parameters chosen 
for our gedanken experiment. Our solution is valid if rj >> (Df3) 2 / V8. Set- 
ting T = T max or D(3rfl 2 = 1, we obtain if >> 1/VS, which is definitely 
true, assuming rj » 1. Next, the condition of the validity of the high tem- 
perature approximation is D » 1. For T = T max , it follows that if'/ 2 « Q. 
This inequality is roughly satisfied [rf/ 2 = 1700, Q = 6700). Finally, as we 
mentioned in the Introduction, the master equation fails at times t < h/k B T. 
Thus, the time considered, r = 2 1 / 4 / v /r/, must be much greater than 1/D, 
which is definitely wrong. Thus, our condition (54) for the creation of MSCS 
is not justified for the parameters considered. 

Next, we discuss at what values of the parameters a MSCS can be gen- 
erated. First, we emphasize the qualitative difference between two condi- 
tions. 1) The condition for distinguishing two positions of the cantilever. 
2) The condition for distinguishing two position of the cantilever and the 
coherence between these two positions. The first condition is relatively 
simple: T < T max = F 2 /k B k c for equilibrium positions (r >> Q) and 
T < T max = 4QF 2 /irk B k c for r = n. The obvious way to increase T max 
is by decreasing the spring constant k c or increasing the magneto-static 
force F. For r = 7r, one additional way is to increase the quality factor 
Q. Condition 2) for generating an MSCS can be satisfied at temperature 
T < T max = h 7/4 Qkl 3/8 /k B F 3/2 m 7/8 , where m is the effective mass of the 
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cantilever, m = k c /u 2 . At T = T max , the MSCS will be generated if we 
satisfy the inequalities 1 <C rf <C Q, or 

mF 2 

1 « « Q. (55) 

One can see that the regime considered in our paper does not allow free 
manipulation of any parameter but Q. Increasing Q, we can satisfy the 
second inequality and, at the same time, increase T max . In our gedanken 
experiment, the tenfold increase of Q (Q — 67000) provides the validity of 
the right-hand inequality in (55) and increases T max to 3 \xK. 



6 Conclusion 

We have shown that the maximum temperature for a single spin measure- 
ment in MFM can be increased by a factor of Q if one utilizes the initial 
transient process instead of the static displacement of the cantilever tip. We 
have obtained an exact analytical solution of the master equation, which de- 
scribes the Q-times magnification of the maximum temperature. In addition, 
we have found the conditions for generation of macroscopic Schrodinger cat 
state in MFM. 
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Figure 1: Geometry of the proposed gedanken experiment. m F and m p are 
the magnetic moments of the ferromagnetic particle and paramagnetic atom. 
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Figure 2: Schematic view of the Gaussians representing the diagonal elements 
|/0-i/2,-i/2| an d |/0i/2,i/2| (seen from the top) in the (-2,-2') plane. We show 

the centers M and M ++ , the variances a' d (transverse) and a d (parallel), 

and the distance between diagonal centers A<j. 
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